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We discuss the equilibrium properties of dilute Bose gases using a non-perturbative formalism 
based on auxiliary fields related to the normal and anomalous densities. We show analytically that 
for a dilute Bose gas of weakly-interacting particles at zero temperature, the leading-order auxiliary 
field (LOAF) approximation leads to well-known analytical results. Close to the critical point the 
LOAF predictions are the same as those obtained using an effective field theory in the large-A^ 
approximation. We also report analytical approximations for the LOAF results in the unitarity 
limit, which compare favorably with our numerical results. LOAF predicts that the equation of 
state for the Bose gas in the unitarity limit is E/(pV) = 1, unlike the case of the Fermi gas when 
E/{pV) = 3/2. 
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I. INTRODUCTION 

Recently we introduced a new theoretical framework 
for the study of a dilute gas of Bose particles with 
tunable interactions [1] based on a loop expansion of 
the one-particle irreducible (1-PI) effective action in 
terms of composite-field propagators. The auxiliary field 
(AF) formalism makes use of the Hubbard-Stratonovitch 
transformation [2] to rewrite the Lagrangian in terms of 
auxiliary fields related to the normal and anomalous 
densities. Employing general quantum field theoretical 
methods [3-7], the AF formalism is part of a continuing 
effort in the community to apply methods traditionally 
used in high-energy physics[8-18] to the study of ultra- 
cold atomic gases[19]. 

For an interacting dilute Bose gas, the leading- 
order auxiliary field (LOAF) approximation is a non- 
perturbative, conserving and gapless approximation that 
describes a large interval of values of the coupling 
constant, satisfies Goldstone's theorem and yields a 
second-order phase transition to a Bose-Einstein conden- 
sate (BEC) regime. In contrast with other resumma- 
tion schemes, such as the large-A^ expansion[20, 21] or 
the functional renormalization techiques[16-18], here we 
treat the normal and anomalous densities on equal foot- 
ing. LOAF produces the same slope of the linear depar- 
ture of the critical temperature from the noninteracting 
limit derived by Baym et aZ.[21] using a large-N expan- 
sion for the critical theory. Unlike the large-N expansions 
developed by Baym et al., the LOAF approximation can 
be used at all temperatures. Furthermore, one can sys- 
tematically improve upon the LOAF approximation by 
calculating the 1-PI action order-by-order corrections. 
The broken U{\) symmetry Ward identities guarantee 
the preservation Goldstone's theorem order- by-order [5]. 

The detailed derivation of the LOAF approximation 
was discussed recently for the case of dilute Bose [22] and 
Fermi[23] atomic gases. Unlike the case of Fermi gases 
where the LOAF approximation is equivalent to the stan- 



dard Bardeen-Cooper-Schrieffer (BCS) ansatz [8, 9], in 
the case of Bose gases the LOAF approximations leads to 
yet unexplored possibilities. Therefore it is important to 
study analytically the LOAF predictions for Bose gases in 
limiting cases such as the case of weakly-interacting sys- 
tems and also in the unitarity limit, which corresponds 
to the strongly-interacting regime where the s-wave scat- 
tering length, oo, is much larger than the inter-particle 
distance. In the unitarity limit the properties of the sys- 
tem have a universal character [24]. The intrinsic non- 
perturbative character of the LOAF approximation may 
become particularly relevant because the development 
of novel cold atom technology that produce stable, flat 
potentials bound by a sharp edge[25, 26] leads to the 
prospect of studying finite temperature properties of di- 
lute gases, such as the BEC transition temperature, T^, 
superfluid to normal fluid ratio, depletion, and specific 
heat, at fixed particle density p. 

In this paper we focus on the study of the LOAF pre- 
dictions in the broken-symmetry phase. We make contact 
with existing analytical approximations in the weakly- 
interacting limit, such as those discussed in the textbook 
of Fetter and Walecka[27] and the analytical results ob- 
tained close to the critical temperature by Kita[28-30] 
using the related Luttinger-Ward functional. We will also 
show that the analytical techniques developed here can 
be used to study the LOAF predictions in the unitarity 
limit. 

This paper is organized as follows: In Sec. II, we briefly 
review the derivation of the LOAF equations. The LOAF 
effective potential and the derivation of thermodynamic 
properties are outlined in Sec. III. In Sec. IV we special- 
ize to the study of the interacting Bose properties in the 
broken-symmetry phase. In Sec. V we study the zero- 
temperature LOAF results in the weakly-coupling limit 
and we compare with the weakly interacting Bose gas 
theory discussed by Fetter and Walecka[27]. The zero- 
temperature analysis suggests the scaling of the LOAF 
equations discussed in Sec. VI. In Sec. VII we discuss the 



LOAF properties close to the critical temperature in the 
weakly-interacting limit and we compare with the results 
obtained by Kita using a related approximation [30]. An- 
alytical approximations of the LOAF predictions in the 
unitarity limit are discussed in Sec. VIII. We conclude in 
Sec. IX. 



II. LEADING ORDER AUXILIARY FIELD 
(LOAF) FORMALISM 

The detailed derivation of the LOAF approximation 
for the case of dilute Bose gases was discussed recently 
in Ref. 22. For completeness, we will review next the 
salient aspects of the AF-formalism derivation. 

In dilute bosonic gas systems, the classical action is 
given by 
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with dx = dtd'^x and the Lagrangian density 
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Here, fj, is the chemical potential and Aq is the bare cou- 
pling constant. In the auxiliary field formalism we use 
the Hubbard-Stratonovitch transformation [2] to elimi- 
nate the quartic interaction in Eq. (2) by introducing the 
real and complex auxiliary fields (AF), x(^) ^^^ ^(2;), 
related to the normal and anomalous densities. We add 
to Eq. (2) the AF Lagrangian density[3-5] 
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Then, the action becomes 
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dx dx' (l)a{x) G-i"t[x, A]{x, x') (t)''{x') 
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where we introduced the notations 7 — fi^ /{2m) and 

x'^V2x{x)-fi, (6) 



together with a two-component notation, (j)'^{x) = 
{ (l){x) , (j)* (x) } , for a = 1,2. $(a;) and J{x) signify the 
five-component fields and currents. 

The generating functional for connected graphs is 
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with 5[$; J] given by Eq. (4). Performing the path inte- 
gral over the fields (pa, we obtain the effective action for 

X,A,A*, as 
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where 



xdx)^{xix),Aix)/V2,A*ix)/V2}, (9) 

SHx) = {s{x),S{x)lyj2,S*[x)ly/2) . (10) 

Next, we expand the effective action about the station- 
ary points, Xo(a;), defined by bS^?i\x\3\I^Xt{x) = 0. We 
obtain the "gap" equations: 
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where we introduced the notations (/)q[xo](2^) as 

-^okoK^:) = I dx'G[xo]\{x,x')j\x') . (13) 



Both xq and Aq include self-consistent fluctuations. 

Expanding the effective action about the stationary 
point, we write 



Scs[x\J] = Scfi[xQ;J] + 
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where D~, [x, x') is given by the second-order derivatives. 
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evaluated at the stationary points. By keeping the gaus- 
sian fluctuations and Legendre transforming, the one- 
particle irreducible (1-PI) graphs generating functional 
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is the negative of the classical action plus self-consistent 
one-loop corrections in the 4>a and Xi propagators. The 
last term in Eq. (16) is next-to-leading order [5], and is 
not included in the leading-order auxiliary field (LOAF) 
approximation. Hence, the static part of the effective 
action per unit volume is 
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In the imaginary time formalism, the last term in Eq. (17) 
becomes 



^Trln[G-Mx]]=/7^|'^ + -m 
where the dispersion relation is given by 
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with ek = 7^^ . At the minimum of the effective potential, 
we also have 
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Using the U{1) gauge symmetry, we choose 0o to be real 
in the broken-symmetry phase. Then, A is real and the 
dispersion, w^ = ^ki^k + 2A), represents the Goldstone 
theorem. 



III. THERMODYNAMICS IN THE LOAF 
APPROXIMATION 

Using standard regularization techniques for effective- 
field theories [31, 32], the renormalized effective potential 
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represents the grand potential per unit volume, 

V,ff = n[T,t^,V]/V. (22) 

Here, po — 4'o i^ the condensate density, and the renor- 
malized coupling constant is related to the s-wave scat- 
tering length by A = 8717 oq . The values of x' and A are 
obtained by solving self-consistently the gap equations 



P = Pq + 



A 



Po = 



(2^)3 


/Cfc + X' 


d^k 


(^ 



{2^Y\2uk 



1 
" 2 

2€k) 



dPk efc + x' 
(27r)3 



^k 



(23) 



d^k A 

(2^^ 



n^, (24) 



where n^ = 1/(6"*=/^ — 1) is the Bose-Einstein particle 
distribution and we assume fcs = 1 units. 

Using the grand potential, 0[r, /^, V], we calculate the 
total number of particles 



the pressure 

entropy 

The energy is obtained as 

E = n + TS + ^iN . 
Hence, the physical density is given by 

P = -9^Kff = ^ (>^' + ^) ' 

the pressure is 

P = -Vcs[po,p.,T] , 

and the entropy density, s = S/V, is 
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The energy density, e = E/V, is obtained as 

e = -p + T.s + iip. (32) 

IV. BROKEN SYMMETRY PHASE 

In the following we will focus on the broken symmetry 
region of the phase diagram, i.e. the regime where the 
density of the BEC condensate is nonzero, po 7^ 0. In the 
broken-symmetry phase, T < Tc, we have equal normal 
and anomalous densities, x' = A = A, and the dispersion 
relation becomes 



Wfc = efc(efe + 2A) 



(33) 



We note that in the long-wavelength limit, Eq. (33) re- 
duces to the linear dispersion relation 



cjfc ~ hk 



0, 



(34) 



with the characteristic velocity (zero sound), y/A/m. 
Comparing (33) with Eq. (21.11) in Ref. 27, i.e. 
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we find that the parameter x' = A = A in LOAF plays 
the role of pV{0) in the weakly interacting Bose gas the- 
ory discussed by Fetter and Walecka, with the zeroth 
moment of the potential 
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In the broken-symmetry phase we can calculate explic- 
itly the temperature independent integrals in Eqs. (21), 
(23) and (24), as 
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with S = A/T. Hence, the gap equations (23) and (24) 
become 
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and the pressure (26) reads 
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The entropy density is given by 
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V. ZERO TEMPERATURE PROPERTIES IN 
THE WEAKLY-INTERACTING LIMIT 

With the above results we can study the properties 
of the zero-temperature Bose gas. For T — Q, the gap 
equations combine to give 
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We also have Eq = — po + P-oP^ with /ip = 2Ap — Aq. Here 
we note: Ap = {^.ttCiIP^'^, ^-4 = {2/'n)"f-'^/'^ao, with 
the dimensionless parameter, ^ — p^'^aQ. 

It is convenient to introduce the following rescaled vari- 
ables 
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Iterating Eq. (52), we obtain the following successive 
approximations for Aq, 
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The coefficient of the ^3/2 term above is the same as in 
Eq. (22.20) inRef. 27. Substituting Eqs. (57) in Eq. (53), 




FIG. 1. (Color online) Comparisons of the exact zero- 
temperature values of the auxiliary field, Ao, condensate frac- 
tion, po/p, pressure, po, and their respective first- and second- 
order approximations, as a function of ^ = p^''^ao, in the 
weakly-interacting regime. The LOAF approximations for the 
auxiliary field, Aq, condensate fraction, po/p, and pressure, 
Po, in the weakly- interacting regime, are given in Eqs. (58), 
(61) and (63), respectively. 



Again, the coefficient of the ^^^^ term above is the same 
asinEq. (22.19) in Ref. 27. 

We note that the next-to-leading order correction to 
the zero-temperature energy was calculated by Wu [33], 
yielding a second-order logarithmic term that cannot be 
captured by LOAF, which is only a one-loop approx- 
imation. The second-order logarithmic correction was 
later confirmed by Hugenhltz and Pines [34] and by 
Sawada [35]. 

For illustrative purposes, in Fig. 1 we depict the exact 
zero-temperature values of the auxiliary field, Aq, con- 
densate fraction, p^/ p, pressure, pq, and their respective 
first- and second-order approximations, as a function of 
^ = p^"ao, ni the weakly-interacting regime. 

By construction [22], in the weak coupling limit 
the LOAF approximation agrees with the Bogoliubov 
approximation[36, 37], which represents the leading- 
order low-density approximation of the theory. In ad- 
dition, the related Popov approximation can be obtained 
from Eqs. (23) and (24) by setting A — y^ — Xpo and 
neglecting the quantum fiuctuations in the anomalous 
density. We showed in Ref. 22 that the LOAF and the 
"gapless" Popov approximation results become qualita- 
tively similar in the weak coupling limit, even though the 
order of the phase transitions remains different. As a con- 
sequence, the LOAF results in the weak-coupling limit 
discussed above agree with the Bogoliubov and Popov 
approximations discussed for instance in the Andersen's 
review article [37]. 
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The coefficient of the ^1"^ term above is the same as in 
Eq. (22.14) in Ref. 27. Similarly, from Eqs. (54) and (57), 
we obtain 
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VI. RESCALED EQUATIONS 



The above zero-temperature results suggest that Xp is 
one of the two characteristic energy scales of the BEG 
system. At finite temperature, we find the second en- 
ergy scale is given by Tq = 47r7[p/(^(3/2)]^/'^, the critical 
temperature of the non-interacting Bose gas. It appears 
that To represents the temperature scale of the BEG sys- 
tem, and it is convenient to supplement the set of scaled 
variables given in Eq. (51), by introducing the scaled tem- 
perature, T = T/Tq. In this context, we note the useful 
results Xp/To = c^ and {A/p)To^^ = 2 [V^C(3/2)]~\ 
withc = 2C3(3/2). 

In terms of the rescaled variables, we have 
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&2 = &i • As indicated above, our approximations of the 
integrals (66), (67) and (68) differ from Kita's approxi- 
mations - see Eqs. (48a), (48b), and (48c) in Ref. 30 - 
at the second order in the 62 and 63 coefficients. 
Then, the gap equations (65) read 
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We introduce the rescaled entropy density as 
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VII. CRITICAL PROPERTIES IN THE 
WEAKLY-INTERACTING LIMIT 

In the following we will follow closely the approach 
outlined by Kita in Ref. 30. It is important to note that 
despite the fact that Eqs. (69), (70) and (71) are differ- 
ent from Kita's Eqs. (41)-(43), in the weakly-interacting 
limit they become the same. Therefore, in the broken- 
symmetry phase, for T < Tc, our results match closely 
Kita's results. Differences are simply due to the fact that 
we found better second-order approximations of the in- 
tegrals (66), (67) and (68). 

In the weakly-interacting limit, wc have 5 <^Q and the 
temperature-dependent integrals can be approximated as 
(see App. A) 
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with c = 2(3(3/2). Combining Eqs. (79) and (80), we 
obtain 
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At the critical point we have po = 0. Then, from Eq. (82), 
we obtain 
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The expansion of T^ in powers of ^ is obtained from 
Eqs. (81) and (83) as 
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Here, the linear coefficient is « 2.33, whereas the 
quadratic coefficient is w —3.23. The result for the 
linear coefficient in ( is the same as that obtain by 
Baym et al. using the large- iV expansion for the critical 
theory[20, 21]. Kita also obtained this linear coefficient, 
see Eq. (52) in Ref. 30. 

We note that the large- A^ expansion results obtained 
by Baym et al. for the critical theory[20, 21] were later 
improved by Kleinert[38] and Kastening[39, 40] using 
five-, six- and seven-loop variational perturbation the- 
ory, respectively. At the seven-loop order, Kastening[40] 
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FIG. 2. (Color online) Comparisons of the exact critical values 
of the auxiliary field, Ac, the ratio {Tc — To)/To, and pressure, 
Pc, and their respective first- and second-order approxima- 
tions, as a function of ^ = p^'^ao, in the weakly-interacting 
regime. The LOAF approximations for the critical values of 
the auxiliary field, Ac, the ratio (Tc — To)/To, and pressure, 
Pc in the weakly-interacting regime, are given in Eqs. (87), 
(85) and (92), respectively. 



calculated a value (T^ - To)/To = (1.27 ± 0.11)^, which 
is in excellent agreement with Monte Carlo lattice field- 
theory results[41-43]. The quadratic coefficient was cal- 
culated by Arnold, Moore and Tomasik[44] , yielding also 
a second-order logarithmic correction that cannot be cap- 
tured by LOAF, which is only a one- loop approxima- 
tion. Their result, (T^ - To)/To « (1.32 ± 0.02)^ + 
[19.7518 ln(0 + (75.7 ± 0.4)]^^^ indicates that two-loop 
contributions are also important in determining the value 
of the quadratic coefficient, as the LOAF result is too 
small. 

Substituting Tc from Eq. (84) in Eq. (83), we obtain 



= hc^ + bic{lblc-b'^c]C 
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In the latter, the quadratic coefficient is 87r/C3(3/2) w 
13.25, whereas the coefficient of the cubic term is w 
—34.76. The leading-order approximation was also ob- 
tained by Kita, see Eq. (53) in Ref. 30. 

The temperature dependence of A close to Tc is derived 



xiT) - x{Tc) 



= -b,b2cH+\bjcH^ + -c4(l-^ 



T 



Tc 



as 



A(r) A(T,) 



To 



To 



= blc^ C' + \/ 46lc4 ^4 + 662c3 e^ 1 



(89) 



T 
Tc 



+ ■ 



The Eq. (88) above is the same as Kita's Eq. 54. Simi- 
larly, we obtain 



P 



^blc^+Jb\cH^ + -blc£,{l- 



T 
T~c 



(90) 



The critical pressure in the weakly-interacting limit is 
obtained from 



Pc/p ~ 1 



This gives 



T 



5/2 



2^Cn3/2) 



{b'o 



46i 



36i 4 



Pc/P = 



C(5/2) 1 
2Ct(3/2)C 
blc" 



5 C(5/2) 



bjc 



\ 5 C(5/2) 
2Ci(3/2) I 3C(3/2) 



6cS(3/2) 

l-ih + Kj]]^- 



(91) 



(92) 
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In the noninteracting limit, we obtain Pc{To)/{jp^'^) = 
47rC(5/2)/Ct(3/2) « 3.4. The coefficient of the linear 
term is 87r[l -|- (107r/3) C(5/2)/C3(3/2)] w 44.94, whereas 
the coefficient of the quadratic contribution is « —6.42. 
For illustrative purposes, in Fig. 2 we depict the exact 
critical values of the auxiliary field, Ac, the ratio {Tc — 
Tq)/Tq, and pressure, pc, and their respective first- and 
second-order approximations, as a function of ^ = p^'^ao, 
in the weakly-interacting regime. 



VIII. CRITICAL PROPERTIES IN THE 
UNITARITY LIMIT 

The LOAF approximation is a non-perturbative ap- 
proximation. Therefore one can solve Eqs. (69) and (70) 
for arbitrary values of coupling constant related to the 
dimensionless parameter, ^ = p^'^ao, to calculate the 
values of the condensate fraction, auxiliary field. A, and 
all thermodynamical variables derivable from the grand 
potential related to the pressure (71). For illustrative 
purposes, in Fig. 3 we depict the zero-temperature val- 
ues of the auxiliary field Aq, condensate fraction, po/p, 
and chemical potential ^O) as a function of ^, whereas in 
Fig. 4 we show the dependence of the critical values of the 
auxiliary field, Ac, the ratio (Tc — To)/To, and chemical 




10^ 103 




FIG. 3. (Color online) Zero-temperature values of the auxil- 
iary field, Ao, condensate fraction, po/p, and chemical poten- 
tial, /io, as a function of ^ = p^'^ag. The condensate fraction, 
po/p, and chemical potential, po/{\p), plots depicted in the 
bottom two panels, were first shown in Fig. 2 of Ref. 22. The 
reader is directed to Ref. 22 for further studies and deriva- 
tions. Here, we emphasize that the numerical values of the 
auxiliary field, Ao, condensate fraction, po/p, and chemical 
potential /^o, in the unitarity limit, compare well with the ex- 
act solutions given in Eq. (96), (97) and Eq. (98), respectively. 



FIG. 4. (Color online) Critical values of the auxiliary field, 
Ac, the ratio {Tc — To) /To, and chemical potential, /ic, as a 
function of ^ = p^'^ao- The critical auxiliary field, Ac, and 
the ratio {Tc — To)/To, depicted in the top two panels, were 
reported previously in Fig. 6 of Ref. 22, and further discus- 
sions regarding these plots can be found in Ref. 22. Here, we 
emphasize that the numerical values of the critical auxiliary 
field, Ac, and the ratio {Tc ~ To)/To, in the unitarity limit, 
compare well with the semi-analytical approximations given 
in Table 1. The numerical value of the chemical potential ^Lc, 
agrees with the exact result obtained in Eq. (98). 



potential, /Xc, as a function of f = p^/^aQ. From Fig. 4 
we conclude that in the unitarity limit (i.e. in the limit 
i -^ oo), we have A^/To -^ 0.873 and Tc/Tq -> 1.396. 

In the previous sections we focused on analytic approx- 
imations of the solutions to the gap equations, Eqs. (69) 
and (70), for the weakly-interacting limit. These approxi- 
mations were derived by dropping the terms proportional 
to A3/2 and higher in Eqs. (69), (70) and (71), and by 
approximating the integrals (66), (67) and (68) as de- 
scribed in App. A. This approach is made possible by 
the fact that the critical value of the auxiliary field, A is 
small in the weakly- interacting limit, and A <C 1 leads 
to (5 = A/f < 1 as f > 1. 

A similar approach is possible for obtaining analyti- 
cal approximations for the critical values of the auxiliary 
field, Ac, and the critical temperature Tc, in the unitarity 
limit. This approach requires improving the approxima- 
tions to the integrals (66), (67) and (68) described in 
App. A, by supplementing the approximation of /2(<^) 
by the term (19/64)6i(5'^/^, and by adding the term 
- (99/256) &i (53/2 ^Q ^Yve approximation of Ii{5) - hid). 
This allows us to take into account the terms proportional 
to A3/2 in Eqs. (69), (70) and (71). We wih describe this 
procedure next. 



As already discussed, the unitarity limit corresponds 
to the strongly-interacting limit, S, — >■ c». In this regime, 
the energy scale Xp diverges, and the only scale remain- 
ing in the problem is the temperature scale, Tq. Hence, 
in the unitarity limit the gap equations are scaled by in- 
troducing the variable, A = A/Tq. From Eqs. (44) and 
(45), we obtain 



1 

P 

with S = A/T : 
obtain also 



Pa 
P 



2 [2 A3/2 
3V^C(3/2) 
2 A3/2 



f^/^h{d) 



f^"h(5), 



(93) 
(94) 



IT C(3/2) 
A/T. By adding the gap equations, we 



_ 8 [2 A3/2 



f^/^[h{5) ~ h(5)] 



(95) 



At zero temperature, T = 0, we can solve (95) to ob- 
tain the value of the auxiliary field in the unitarity limit 
(UL) at zero temperature, as 



Ao/To = 



C(3/2) 



2/3 



1.147. 



(96) 
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TABLE I. Exact and approximate values of the critical 
values of the auxiliary field, Ac, and the critical temperature, 
Tc, in the unitarity limit. 

Exact 1=* order 2"^^ order 3'''* order 



VSl 0.790 1.253 
Ac/To 0.873 4.619 
Tc/To 1.396 2.941 



0.694 0.789 
0.653 0.875 
1.356 1.407 



and from (94) we obtain the UL condensate fraction is 



Po 
P 



(97) 



Alternatively, the UL condensate depletion fraction 
is 1/4. 

The critical temperature, 5c in the unitarity limit is 
obtained by solving Eq. (94) with po = 0. Assuming the 
expansion of l2{5) as a series powers of 5 is truncated 
by dropping terms proportional to 5^/"^ and higher, we 
find 5c as the solution of a quadratic equation. Then 
we can calculate fc = Tc/Tq from Eq. (95). Table I 
summarizes approximations for the critical values of the 
auxiliary field, Ac, and the critical temperature Tc, in 
the unitarity limit. These approximations are compared 
with the numerical "exact" values. The errors of the 
third-order approximations relative to the exact values 
are less than half of a percent. 

Finally, we can show that the equation of state in the 
unitarity limit is independent of temperature. We begin, 
by calculating the UL asymptotes of the chemical poten- 
tial and pressure. For arbitrary temperature, Eqs. (56) 
and (71) give 



fi = 2Xp, 



p = Xp' 



(98) 



Then, from Eq. (55) we find the UL equation of state at 
zero temperature is 



Eo/{poV) = eq/po ^ eo/po = 1 



(99) 



For finite temperature we use instead Eqs. (75) and (73) 
and obtain 



E/{pV) = e/p = e/p = 1 



(100) 



This result is different than in the Fermi gas case, 
where the equation of state in the unitarity limit is, 

E/{pV)=3/2 [23, 45]. 



IX. CONCLUSIONS 

In summary, in this paper we discussed analytical ap- 
proximations of the properties of dilute Bose gases us- 
ing the LOAF approximation in the weakly-interacting 
and the strongly-interacting (unitarity) limit. We focus 



deliberately on the case of Bose gases in the broken- 
symmetry phase, in order to make contact with ana- 
lytical results already existing in the literature. The 
weakly-interacting results at zero-temperature are shown 
to be identical with those found in the weakly interact- 
ing Bose gas theory discussed by Fetter and Walecka[27], 
whereas close to the critical temperature, the LOAF re- 
sults are similar to those obtained by Kita using the re- 
lated Luttinger-Ward functional[28-30]. In obtaining our 
results, we have improved the analytical approximation 
of temperature-dependent integrals, first discussed by 
Kita[30]. These approximations were then applied to the 
analytical study of the LOAF predictions in the unitarity 
limit and found to give good agreement with our numeri- 
cal results. LOAF predicts that the equation of state for 
the Bose gas in the unitarity limit is E/{pV) = 1, unlike 
the case of the Fermi gas when E/{pV) — 3/2. 
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Appendix A: Approximations of certain integrals 

For completeness, in this appendix we will derive 
the first- and second-order approximations of the inte- 
grals (C6), (67) and (68). Our approach follows closely 
the discussion in Kita's paper, see Ref. 30. Our first-order 
approximations are the same as Kita's, but we differ at 
the second order. 

Consider the integral 



hiz) 



e + z x/e 



de 



¥C(3/2)7o w e'^-l' 



(Al) 



with Lj = \/e{e + 2z). In the weakly-coupling limit we 
have z — > 0, as A — > 0, and we seek a power expansion of 
/i(z) in powers of z. We have 



%A _ %A ye 



w e" - 1 e" -1 {e^ - If 



+ 0{e^). (A2) 



Upon integration with respect to e, wc find that only the 
first term converges. 



bn = 



/7?C(3/2) Jo 



Jo e^-1 



(A3) 



Therefore, we can write Ii{z) = bo + Ii_ao{z), with the 
remainder 



h.,ooiz) 



de 



' e + z ^e ^ 



7rC(3/2)yo \ uj e"-l e^ - 1 



(A4) 
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The leading-order approximation of /i,oo(2:) is obtained 
by expanding out the exponentials in the denominator to 
first order. We obtain 



'1,00 



(^) 



7rC(3/2)yo \e + 2z 



\e + 2z J 



de 



'hwfz . 



(A5) 

with h\ = \/27r/C(3/2). Note that it is important to treat 
on equal footing the two exponentials in Eq. (A4) in order 
to obtain the correct divergence subtraction. 

The second-order correction to Ii.oo{z) corresponds to 
the term linear in z in the power expansion of 



z 1 



+ 2z ,/e Vl + i^e(e + 2z) 



This gives 



7rC(3/2)7o (e + 2)- 



;de = b2 z 



(A6) 



(A7) 



with &2 = bi/2. Therefore, the second-order approxima- 
tion of Ii{z) is 



hiz) 



V^ai/2)Jo UJ e 



^e + z ^/e de 



1 



1 - bi ^/z + b2 z . 

(A8) 

The first- and second-order approximations of Ii{z) are 
illustrated in Fig. 5. The coefficient 62 = 0.559 calcu- 
lated by Kita[30] results in a worse second-order approx- 
imations. 

Similarly, we can evaluate the expansion in powers of 
z of the integral 



W) 



/^a3/2)Jo ueJ 



z yfe de 



1 



(A9) 



In this case, there is no analytic part identified after per- 
forming the z expansion of the integrand in l2{z)- There- 
fore, we have l2{z) — l2,oo{z), with the remainder 



Az) 



de 



y^C(3/2) Jo V^T2^ e- 



(AlO) 




"0.6 



FIG. 5. (Color online) Comparisons of the exact 2 depen- 
dence of /i(z), and the first- and second-order approximations 
in powers of z. 
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FIG. 6. (Color online) Comparisons of the exact z depen- 
dence 01/2(2), and the first- and second-order approximations 
in powers of z. 
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FIG. 7. (Color online) Comparisons of the exact z depen- 
dence of 13(2), and the first- and second-order approximations 
in powers of 2. 



The first-order approximation of l2,oa{z) is obtained by 
expanding out the exponential in the denominator to first 
order. We obtain 



^(1) 
'2,00 



(^) 



/^C(3/2) Jo 



de 



2z^ 



b[^, (All) 



with b'^ = bi- The second-order correction to l2,oo{z) cor- 
responds to the term linear in z in the power expansion 
of 



1 



1 



e + 2z V^\l + l^eie + 2z] 



- 1 



(A12) 



This gives 



r(2) 



W = - 



de 



V^C(3/2)7o e + 2yi 



= -b'2Z, (A13) 



with 62 ~ bi. Therefore, the second-order approximation 
oi hiz) is 



hiz) 



/^C(3/2) Jo Lo e' 



z ^e de 



1 



bx^z-\)2Z. (AI4) 



The first- and second-order approximations of 12(2) are 
illustrated in Fig. 6. The coefficient 63 = 1-12 calculated 
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by Kita[30] results in a worse second-order approxima- 
tions. 

Finally, we also need an approximation for the integral 



h{z) 



/^C'(3/2) 



/ Vedeln(l-e-"). (A15) 
Jo 



We note that the derivative of Iziz) with respect to z can 
be written in terms of Ii{z) and l2{z), as 



dhjz) ^ 2 

dz y^C3(3/2) 



(A16) 



Therefore, we can write 



^3(^) = 7^7^ 1^0 + /' rf^' [hiz') ~ h{z')\ } , (A17) 
C-H3/2) '^ Jo J 

where 6q = — C(5/2)/C(3/2). Thus, we obtain the second- 
order approximation of I-i{z) as 



h{z) 



46i 



36i 



3/":)^ V 3 4 



Ct(3/2) 



2"'" + -7^2") • (A18) 



The first- and second-order approximations of I^iz") are 
illustrated in Fig. 7. 
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